Effect of final state interactions on particle production in d+Au collisions at RHIC 
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Abstract 

We show that particle species dependence of enhanced hadron production at intermediate transverse momentum (pr) for d+Au 
collisions at RHIC can be understood in terms of the hadronic rescatterings in the final state. A multiphase transport model 
(AMPT) with two different hadronization mechanisms: string fragmentation or parton coalescence, is used in our study. When 
the hadrons are formed from string fragmentation, the subsequent hadronic rescatterings will result in particle mass dependence 
of nuclear modification factor Rcp, which is consistent with the present experimental data. On the other hand, in the framework 
of parton coalescence, the mass dependence disappears and the strangeness plays an important role. Both mechanisms failed to 



reproduce the pr dependence of Rcp of pion, indicating that initial-state effects might be also important in such collisions. 
, Keywords: Nuclear modification factor, hadronization mechanisms, hadronic rescatterings 
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1. Introduction 

' The Cronin effect, which refers to the enhanced hadron pro- 
duction at intermediate transverse momentum (pr) with in- 
creasing target nucleus size in proton-nucleus (pA) collisions, 
,was first observed by Cronin et al. [1, 2] in 1975 at center 
of mass energy s/snn - 27.4 GeV. Recent experimental data 
in li-i-Au collisions from RHIC have shown that similar ef- 
fect exists at higher collision energy (V^ = 200 GeV) llli. 
An adequate understanding of the Cronin effect becomes espe- 
cially important for making reliable theoretical interpretations 
of the observations in heavy ion collisions at RHIC, in which 
the quark-gluon plasma (QGP) is thought to be created Istl- 
Traditional explanations of the Cronin effect all involve mul- 
tiple scattering of incoming partons that lead to an enhance- 



ment at intermediate pj 1161- 11 ID . The models can reproduce 
the observed centrality dependence for pions very well. How- 
ever, none of these initial-state models would predict a species- 
dependent Cronin effect, as initial state parton scattering pre- 
cedes fragmentation into the different hadronic species |4] . This 
indicates possible final-state interactions (FSI) would possibly 
contribute to the Cronin effect. 

Recently, Hwa and Yang |12] demonstrated the recombi- 
nation of soft and shower partons in the final state could ex- 
plain the mass dependent Cronin effect. This model do pre- 
dict a larger enhancement for protons than for pions at 1 < 
Pt < 4- GeV/c. However, the inclusion of recombination 
from deconfined partons requires a high energy density in ini- 
tial state, which may not be justified in li+Au collisions. An 
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alternate way for hadronization is string fragmentation 1131 - 
1511 . In this scenario, hadrons are formed from the decays of 
excited strings, which result from the recombination of ener- 
getic minijet partons and soft strings that produced from ini- 
tial soft nucleon-nucleon interactions. After the hadronization, 
the followed rescatterings between formed hadrons or between 
formed hadrons and nucleon spectators should also be taken 
into account. 



In this paper we study quantitatively how the two different 
hadronization mechanisms (string fragmentation and parton co- 
alescence) and the followed final state interactions would con- 
tribute to the nuclear modification factors in d+An collisions at 
= 200 GeV. A multiphase transport model (AMPT) lfl4l[l5|] 
with two versions: default (hadronization from Lund string 
fragmentation, mainly hadronic rescatterings in the final state, 
version 1.11) and string melting (hadronization from quark co- 
alescence, version 2.11) are used to study the later-stage ef- 
fect. The final state hadronic and/or partonic interactions are 
included in the calculations. Quark transverse momentum kick 
due to multiple scatterings are treated in the same way as in Ref. 
LI 3]- To have a clean test of final state interactions, we assume 
no quark intrinsic pr broadening due to initial multiple parton 
scattering f^, and see how the final state interactions would 
contribute to the Cronin effect observed. We show that recent 
data on particle species dependence of central-to-peripheral nu- 
clear modification factor Rcp at midrapidity for tZ-i-Au colli- 
sions at RHIC can be understood in terms of the hadronization 
from string fragmentation and the followed hadronic rescatter- 
ings in the final state. 
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2. The AMPT model 

The AMPT model is a hybrid model that consists of 

four main components: the initial conditions, the partonic in- 
teractions, conversion from the partonic to the hadronic matter, 
and the hadronic interactions. The initial conditions, which in- 
clude the spatial and momentum distributions of hard minijet 
partons and soft string excitations, are obtained from the HI- 
JING model (version 1.383 for this study). One uses a Woods- 
Saxon radial shape for the colliding gold nuclei and introduces a 
parameterized nuclear shadowing function that depends on the 
impact parameter of the collision. The ratio of quark structure 
function is parameterized as the following imp act-p arameter- 
dependent but (and flavor)-independent form il4ll 



RA{x,r) = 



1 + 1.191n'/*A(x^ - l.lx^ +Q.2\x) 



1.08(Ai/3-l)VI .,0.01 
ln(A +1) 



(1) 



where x is the light-cone momentum fraction of parton a, and 
fa is the parton distribution function. The impact-parameter 
dependence of the nuclear showing effect is controlled by 



aA(r) = 0.133(A 



1/3 , 



1 - r^IRl 



(2) 



with r denoting the transverse distance of an interacting nucleon 
from the center of the nucleus with radius Rq - 1.2A'^-^. The 
structure of deuteron is described by the Hulthen wave function. 
Scatterings among partons are modeled by the Zhang's parton 
cascade (ZPC) [16], which at present includes only two-body 
elastic scatterings with cross sections obtained from the pQCD 
with screening masses. In the default AMPT model, after par- 
tons stop interacting, they recombine with their parent strings, 
which are produced from initial soft nucleon-nucleon interac- 
tions. The resulting strings are converted to hadrons using the 
Lund string fragmentation model. In case of string melting, 
the produced hadrons from string fragmentation, are converted 
instead to their valence quarks and antiquarks. The followed 
partonic interactions are modeled by ZPC. After the partons 
freeze out, they are recombined into hadrons through a quark 
coalescence process. The dynamics of the subsequent hadronic 
matter is described by a hadronic cascade, which is based on a 
relativistic transport model (ART). Final hadronic observables 
including contributions from the strong decays of resonances 
are determined when the hadronic matter freezes out. 



3. Results 



We learn that Lin and Ko have done a nice study ||15|] 
on global properties of deuteron-gold collisions with default 
AMPT model, which shows good agreement with later exper- 
imental data i flTi [Tsll . Their study on nuclear effects is up to 
Pt = 2 GeV/c. Here we focus on the intermediate to higher pr 
range where the Cronin effect exists. We study the deuteron- 
gold collisions at y/s^N - 200 GeV. The string fragmentation 
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Figure 1: Transverse momentum spectra of mid-rapidity (|>'| < 0.5) pions, 
kaons, protons and (p mesons in "minimum bias" d+Au collisions from default 
AMPT (solid lines) and string melting AMPT (dashed lines) versus data from 
STAR Collaboration (statistical error only) 



parameters are chosen to be the same as in Ref. [15]. The par- 
tonic scattering cross section is chosen to be 3 mb. The events 
are separated into different centrality bins using the number of 
participant nucleons suffering inelastic collisions. Fig. 1 shows 
the mid-rapidity {\y\ <0.5) transverse momentum spectra of pi- 
ons, kaons, protons and mesons for "minimum bias" d+Au 
collisions from default AMPT (solid line) and string melting 
(dashed line). It is seen that both default and string melting 
AMPT can reproduce the tt'^ and K"^ spectra well. For proton 
and antiproton production, the default version works well for 
Pt > I GeV/c, but underestimates the low pr proton and an- 
tiproton yields. The string melting version underestimates the 
proton and antiproton production in the whole pr range. For 
(p meson spectrum, the default version works well in the whole 
Pt range, while the string melting one overestimates the low pj 
(p meson yields in the "minimum bias" li-i-Au collisions. 
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Figure 2: (Color online) (a) Rcp for n^, K^, p, (p, A and H in d+Au yijviv 
= 200 GeV collisions from default AMPT without final state interactions in- 
cluded, (b) Kaon and A rapidity distribution in "minimum bias" d+Au colli- 
sions from default AMPT without final state interactions included, (c) quark (u, 
u, s, and .v) rapidity distribution before coalescence in "minimum bias" d+Au 
collisions from sfiing melting AMPT. 
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To study the final state effect on the nuclear modification fac- 
tor Rcp, we first calculate the Rcp (0-20%/40-100% centrality) 
of different hadrons without including the final state hadronic 
interactions and resonance decays in the default AMPT. The 
Rcp, which compares particle yield from central collisions to 
that of peripheral collisions, is defined as the ratio of particle 
yields in central collisions over those in peripheral ones scaled 
by number of inelastic binary collisions N^^^, that is, 

[dNI(NbinPTdpT)]cQntm\ 



R 



CP - 



(3) 



[dN/ (NbmPTdpT)]pQj:iphQj:ai 

Here we use the same Nbi„ value as STAR Collaboration at the 
corresponding collision centrality |3]. One can see in Fig. 2(a) 
that there are only slight differences for the Rcp of different 
particle species. This is because hadrons are produced from 
string fragmentation in the Lund model, and the fragmenta- 
tion patterns for different particle species in central collisions 
and in peripheral collisions are set to be the same. We note 
that the Rcp of proton and A would be larger due to the asso- 
ciate production A^-i-A^ N+A+K'^ from initial-state multiple 
scatterings in the gold beam direction (y < 0) at large rapid- 
ity, as shown in Fig. 2(b) the enhanced production of A and 
K'^ in "minimum bias" d+Au collisions. Some of the parti- 
cles are scattered into mid-rapidity region. The strange quark 
enhancement in the large rapidity region (gold beam direction) 
will cause the corresponding increase of s quark at other rapid- 
ity regions due to net strangeness conservation. This is shown 
in Fig. 2(c): the quark rapidity distribution before coalescence 
in "minimum bias" d+An collisions from the string melting 
AMPT. 
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Figure 3: Mid-rapidity Rcp Qy\ < 0.5) in d+Au ^fsNN = 200 GeV collisions 
from default AMPT witii (red circles) and without (black circles) final state 
interactions. 



and resonance decays. For jt^ and K , the Rcp decreases after 
including the final state interactions from intermediate to high 
Pt- And this suppression increases with pr- Since most of the 
produced particles are pions at mid-rapidity, the scatterings are 
mainly the particle-pion interactions for p » mo, where niQ is 
the mass of corresponding particles. For n-n elastic collisions, 
the resonance peak centers at the position of n - n center-of- 
mass energy yfs^ close to p meson rest mass. Since most of 
the outgoing particles which probably scatter with each other 
are in the similar directions, the open angle between two scat- 
tering particles is small. In our studied pr and rapidity range, 
for one particle at low pr, and another particle with higher en- 
ergy, the calculated -\fs^ is closer to the resonance peak. As a 
result, this hadronic effect on Rcp is enhanced with pr for pion. 
Similar argument is also valid for K-tt scatterings. For heavier 
particles like antiproton, (f> meson. A, the Rcp increases at low 
Pt according to their corresponding production channels or due 
to diffusions into mid-rapidity region, but changes slightly for 
Pt > 3 GeV/c. For even heavier particle, like S , there is no 
obvious change of Rcp due to the final state interactions. Here 
particle mass plays important roles in the hadronic rescatter- 
ings. It will determine the space-time configuration of formed 
hadrons 02011. 
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After including the final-state hadronic interactions and 
strong decays of resonances, the Rcp of different particle 
species will change differently, as they have different masses 
and scattering cross sections. We show in Fig. 3 the com- 
parisons of Rcp with and without the final state rescatterings 



Fi gure 4: (Color online) Mid-rapidity (|y| < 0.5) Rcp iti d+A^u ^ — 
200 GeV collisions: (a) default AMPT with final state interactions; (b) 
string melting AMPT with final state interactions; (c) experimental data of 
Rcp from STAR Collaboration (statistical error only) (J. (d) The ratios of 
Rcp(p)IRcp(J^') from STAR Collaboration, from default AMPT, from string 
melting AMPT, and from default AMPT without final state interactions. 
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According to above analysis, the final state hadronic rescat- 
terings will lead to particle mass dependence of Rcp- In Fig. 4, 
we compare the data with model calculations. At intermediate 
Pt, the Rcp of heavier particles like antiproton, (fi meson. A, 
S will be larger than those of n" and .The result is quali- 
tatively consistent with experimental data |3], as shown in Fig. 
4(c). At intermediate pj, the Rqp of antiproton is systematically 
larger than that of . In Fig. 4(d), the ratio Rcpip)/Rcpi^ ) 
from the default AMPT model also agrees very well with ex- 
perimental data. We note that the present calculation without 
including the possible quark intrinsic pr broadening at initial 
state can not reproduce the pr dependence of Rcp- This dif- 
ference between data and model results may indicate that the 
initial state effect is important. The year 2008 data of RHIC 
with higher statistics will provide more precise measurements 
and test our predictions for other hadron species like (f> meson, 
A, S , etc. 

For comparisons, the Rcp from string melting AMPT with 
quark coalescence is also studied. We have shown in Fig. 
2(c) that the excess of s quark over s quark at mid-rapidity 
is partly due to associate production from initial multiple in- 
teractions. Combining this effect with the coalescence of par- 
tons, there are enhancements of corresponding hadrons at in- 
termediate Pt- The Rcp values for different particle species 
that contain different number of S-quarks are shown in Fig. 
4(b). Note, multi-strange hadrons are particularly interesting 
as they suffer much less hadronic interactions 112 ill compared 
with non-strange hadrons. Therefore they are more sensitive 
to early stage dynamics. At intermediate pr, there is an en- 
hancement of Rcp according to the number of s quarks, that is, 
Rcp(A) < RcpC^^) < Rcp(^ )■ Note that the values of Rcp 
for strange particles (A, H, and Q) are close to each other (not 
shown here) at the same transverse momentum region. If one 
assumes the validity of the coalescence approach, this observa- 
tion shows that the measured Rcp can, to some extend, reflect 
the density of quarks shortly before the freeze-out. However, 
from the coalescence calculation the Rcp of antiproton is close 
to that of TT" at intermediate pr, which is not consistent with 
experimental data, as shown in Fig. 4(d). 

4. Summary 

In summary, we studied the mechanism of hadron formation 
and subsequent interactions in d+Au collisions at = 
200 GeV. In a multiphase transport model with Lund string 
fragmentation for hadronization and the subsequent hadronic 
rescatterings included, we find particle mass dependence of 
central-to-peripheral nuclear modification factor Rcp- Recent 
data on particle species dependence of Rcp at mid-rapidity for 
d+Au collisions at RHIC can be understood in terms of this 
final-state hadronic rescatterings. This shows the importance of 
final state hadronic interactions in cZ-i-Au collisions, since none 
of the initial-state models would predict a species-dependent 
Rcp at present. However, the calculations can not reproduce 
the Pt dependence of Rcp with only final state interactions, 
this indicates the initial state effects might be also important. 



Issues associated with the initial condition such as gluon 
saturation |22-24], parton intrinsic pp broadening |6-ll]] 
and so on were not addressed in this paper. The future more 
complete analysis should include these effects. In comparison, 
if the hadron is formed from quark coalescence, it is difficult 
to explain antiproton transverse momentum spectra and the 
particle species dependence of Rcp- On the other hand, the 
strangeness effect plays an important role. More precision data 
in the future wiU test the findings in this paper. 
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